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rhe  following  report  covers  the  period  from  1  August  1972  to  15  March  1973. 
During  this  period  effort  has  been  concentrated  in  the  following  areas: 

I.  Crystal  Growth 

II.  Magnetic  Resonance  Study  of  Defects 

III.  Hydrogen  Diffusion  Experimental 

IV.  Hydrogen  Diffusion  Theoretical 

V.  Determination  of  Conduction  Band  Parameters  by  Use  of  Hydrogen 
Solubility  Theory  Developed  in  Previous  Report 

VI.  Far  I.R.  Spectrometer 

The  four  preprints  included  in  ‘•he  last  progress  report  have  all  been 
accepted  for  publication  and  will  appear  shortly.  In  addition,  three  papers 
dealing  with  work  in  progress  were  presented  at  the  San  Diego  APS  Meeting. 


I.  Crystal  Growth 

The  relation  between  the  actual  temperature  profile  inside  the  chemical 
vapor  transport  ampoule  and  the  temperature  distribution  of  the  exterior  heat 
source  has  been  followed  further.  The  temperature  gradient  at  the  growing 
crystal  determines  to  a  large  extent  the  interface  stability.  Hence  this 
investigation  is  basic  to  CVT  crystal  growth  in  general.  Preciouf.  little  had 
been  known  before  our  measurements,  concerning  temperature  oscillations  and 
profiles  in  this  crystal  growth  technique.  By  expanding  the  study  from  TeCl 

4 

to  iodine,  that  is  also  frequently  ased  as  CVT  transport  agent,  we  were  able 
to  verify  chemical  convective  instability  on  a  broader  basis,  leading  to  a 
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publication  that  will  be  basic  for  future  CVT  studies.  (See  the  attached 
Preprint  I  that  has  been  submitted  to  the  Journal  of  Crystal  Growth;  an  invited 

talk  on  this  subject  will  be  presented  at  the  Gordon  Research  Conference  on 
Crystal  Growth  in  July.) 

The  crystal  growth  efforts  are  summarized  in  the  er.'los«d  thesis  (Preprint 
II)  of  M.  C.  DeLong.  Beyond  this  thesis,  further  insight  has  been  obtained 
into  the  actual  transport  chemistry  of  the  process.  In  former  progress  reports 

it  had  been  pointed  out  already  that  the  transport  reaction  assumed  by  earlier 
workers 


Ti02(s)  +  TeC14(g)  +  TiC14(g)  4  Te02 (g)  (1) 

is  at  best  a  crude  overall  equation  for  a  chain  of  subreactions.  One  of  these 
reactions,  thermal  dissociation  of  TeCl4  has  been  identified  earlier  in  this 
work.  Quantitative  measurements  of  the  dissociation  kinetics  arm  now  imder 
way.  In  addition,  thermo-chemical  calculations  predict  that  Te02  should  also 
dissociate  to  a  large  extent  in  the  employed  temperature  range.  This  species 
will  be  included  into  the  kinetic  studies  which  will  result  in  better  defined 
thermo -chemical  d-»*a  for  the  various  tellurium  compounds  involved,  than  are 
currently  available.  Furthermore  we  found  that  Si02  of  the  quartz  container 
is  also  transported  in  the  system,  leading  under  certain  conditions  to  a 
sizeable  contamination  of  the  Ti02  crystals.  We  assign  the  observed  Si02 
transport  preliminarily  to  the  coupled  reactions  [that  are  also  coupled  to  (1)]. 


Si°2(S)  *  TlC14(g)  »  Si“4(s1  *  ™2(5) 
SiC14(g)  *  °2(g)  *  Sl02(s)  *  2C12(g) 


(2) 

(3) 
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Further  experiments  arc  planned  to  verify  the  role  of  (2)  and  (3)  which  arc 
of  general  interest  since  various  other  workers  employ  quartz  ampoules  in 
CVT  processes  that  involve  SiCl,. 

4 

II.  Magnetic  Resonance 

Work  in  this  area  has  ccnterd  on  investigation  and  identification  of 
various  charge  states  of  impurities  in  rutile.  See  attached  preprints 
III  and  IV.  In  addition  to  the  work  discussed  in  the  preprints  we  arc 
studying  W-doped  and  V-doped  material  and  are  attempting  to  identify  the 
site  of  interstitial  protons  in  tho  rutile  lattice. 

III.  Hydrogen  Diffusion  experimental 

A  large  effort  has  been  made  in  this  area  because  of  the  importance 
of  H  as  a  probe  in  the  study  of  electron  trapping  levels  as  discussed  in  IV. 
This  has  proved  to  be  an  extrer-ly  difficult  project  both  experimentally  and 
theoretically  (Sec.  IV).  In  order  to  avoid  large  electric  field  effects  :',t  is 
necessary  to  perfor-  the  experiments  at  nearly  constant  Fermi  Level.  This  has 
been  achieved  by  simultaneously  diffusing  Hin  and  D out  or  vice  versa  in  su:h  a 
way  that  the  total  concentration  of  II  ♦  D  remains  nearly  constant.  Problems 
of  vacuum  contamination  at  the  elevated  temperatures  required  have  been  severe 
as  have  surface  effects.  Thus  far  satisfactory  results  have  not  been  obtained 
but  we  feel  that  the  problems  are  nearly  under  control. 
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IV.  Hydrogen  Diffusion  Theoretics! 


A.  Introduction 

Diffusion  in  insulators  is  generally  accom>anied  by  changes  in  the  Fermi 
Level,  which  in  turn  results  in  electric  fields  which  influence  the  diffusion 
rate  in  the  material.  Since  the  electric  field  present  may  vary  greatly  with 
the  sample  and  the  experimental  conditions  one  cannot  expect  a  unique  diffusion 
coefficient  or  even  that  the  system  will  obey  Fick*s  Law.  The  observed  diffusion 
will  be  the  result  of  many  factors  including  microscopic  jump  rates  in  the 
absence  of  fields,  mpuritics,  and  nearby  diffusing  ions  etc.;  distribution  of 
impurities  and  defects;  in  addition  to  internal  electric  fields.  Yet  experi¬ 
mentally  what  is  measured  is  simply  concentration  vs.  time  -  indeed,  often  only 
total  integrated  concentration  vs  time.  The  purpose  of  this  study  (consisting 
of  the  theory  of  this  paper  and  the  experimental  work  of  Sec.  Ill)  is  to  identify 
precisely  wlat  can  be  determined  from  diffusion  measurements. 

As  a  prototype  system  we  will  consider  the  simultaneous  diffusion  of  H  and 
D  in  rutile  in  such  a  fashion  that  the  total  concentration  of  H  plus  D  remains 
nearly  constant.  As  will  be  shown,  this  permits  a  quantitative  treatment  of 
the  effects  of  both  concentration  gradients  and  electric  fields  in  a  self- 
consistent  manner,  and  determination  of  an  effective  diffusion  coefficient 
in  the  absence  of  a  field. 

In  summary,  the  goal  of  this  study  if  to  indicate  the  conditions  under  which 
Fick's  Law  is  valid  and  to  determine  the  meaning  of  the  measured  diffusion 
coefficient  if  these  conditions  are  met. 
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•  Diffusion  Equation 


We  begin  by  deriving  an  equation  describing  diffusion  in  the  presence  of 
both  a  concentration  gradient  and  electric  fields,  in  a  form  suitable  to  the 
present  problem  as  indicated  in  the  Introduction.  We  shall  not  be  concerned 
with  a  first  principles  calculation  of  diffusion  rates  and  in  fact  we  wish  to 
formulate  the  problem  is  a  manner  which  is  as  model  independent  as  possible. 

We  assume  that  the  concentration  n  of  a  given  species  can  be  represented  as 
a  continuous  function  of  position  which  is  sufficiently  slowly  varying  that 
An/n  «  1  (where  An  is  the  change  in  n)  over  n  region  large  enough  that  a 
macroscopic  electrostatic  potential  can  be  used,  i.e..  over  a  region  of  several 

5-10)  atomic  planes.  We  suppose  further  that  the  macroscopic  potential,  $,  is 
sufficiently  slowly  varying  over  the  same  region  that  qA<J>  «  kT  (where  q  is  the  charge 
of  the  diffusing  species  and  A*  is  the  change  in  <J>).  and  that  the  local  concentration 
influences  the  average  (over  the  above  region)  probability  that  a  particular  ion 
will  jump  out  of  the  region  only  through  the  macroscopic  potential  gradient. 

The  latter  is  an  assumption  that  the  concentration  is  low  enough  that  "most" 
hops  of  "most"  ions  are  not  influenced  by  the  presence  of  others .  The  fact  that 
they  may  be  influenced  occasionally  is  covered  by  the  averaging  process,  but 
we  are  ignoring  any  dependence  of  the  average  jump  frequency  on  concentration. 
Considering  the  normal  concentration  range  of  ^  1018-ip.19/Cm3  this  should  be  a 
good  approximation. 

Using  these  assumptions,  and  considering  one-dimensional  diffusion,  it 
is  readily  seen  that  the  time  rate  of  change  of  concentration  is  given  by: 
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3n(x,t) 

"  V  ■  "<x  *  Ax.t)  fL(*  *  Ax)  *  n(x  -  Ax,t)  F.(x  -  4x) 


-  n(x,t)[FL(x)  .  Fr (x) ]  , 

where  x  Is  the  position,  t  is  the  time,  FL  end  FR  are  the  average  probabilities 
of  an  ion  leaving  the  region  by  hopping  to  the  left  and  right  respectively. 

The  relation  between  F^  and  FR  can  be  found  by  considering  a  concentration 
profile  corresponding  to  equilibrium  in  the  presence  of  a  macroscopic  potential, 

♦  Cx).  (Note  that  equilibrium  requires  a  uniform  chemical  potential,  not  a 
uniform  concentration.)  Since  the  F's  depend  on  n  only  through  *  (as  assumed 
above)  the  relation  found  at  equilibrium  will  be  valid  for  all  n.  The  equilibrium 
concentration  is  given  by  Boltzmann  statistics  as: 

n(x)  *  n0e'^W/kT  ,  (2) 

where  nQ  is  a  constant,  k  is  Boltzmann's  constant,  and  T  is  the  absolute 
temperature. 


Since  in  equilibrium  the  net  number  crossing  any  plane  must  be 


zero  we 


have: 


n  e-q4>(x+Ax)/kT  «  r  .Am  -q<f>(x)/kT 
V  FL(x  *  to)  *  n0e^vw,K1  Fr(x)  . 


Expanding  the  exponents  and  rearranging,  we  find: 


FL(X  .  Ax)  %  Fr(x)  [l  .  £  4xJ  . 
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Substituting  (3)  in  (1)  we  find,  after 


some  manipulation,  the  equation: 


where  Ax  is  the  width  of  the  region.  Ax  takes  the  place  of  the  lattice  constant 
in  the  normal  derivation  of  Fick’s  Law.  In  our  case  Ax  is  not  the  lattice 
constant  and  FR(x)  is  not  the  probability  of  jumping  from  one  plane  to  the 
next.  We  can,  however,  combine  them  into  an  effective  diffusion  constant  as 
shown  below.  For  the  momei.t  we  simply  let 

(Axj2  FrW  *  «(*)•  (5) 

What  we  have  done  is  consider  a  region  large  enough  that  local  fluctuations 
(e.g.,  is  there  another  impurity  or  lattice  defect  near  the  diffusing  ion?)  average 
out.  Hence  Ax  does  not  correspond  to  a  single  jump  but  rather  to  many  jumps. 
Exactly  how  many  is  not  important  so  long  as  it  is  enough  to  eliminate  local 
perturbations.  Similarly  FR  is  not  the  usual  jump  frequency  but  rather  1/t  where 
t  is  the  average  time  required  for  a  particular  ion  to  leave  the  region,  Ax,  by 
jumping  to  the  right.  In  the  ideal  case  of  a  completely  homogeneous  crystal, 

with  low  concentration  of  diffusing  ions,  and  in  which  only  single  jumps  occurred, 
we  weald  have: 


where  a  is  the  lattice  constant,  F  is  the  microscopic  jump  frequency,  and  D  is  the 
usual  diffusion  coefficient.  In  general,  however,  none  of  these  conditions  are 
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met  and  g(x)  is  not  related  to  a  and  F  in  any  simple  fashion.  We  will  see 
later  (in  Sec.  D)  the  significance  of  g(x).  Thus: 

-  k0  =  &  [«(x)  *  &  8(x)  n(x,t)  .  (6) 

To  proceed  further  we  suppose  that  g(x)  is  constant.  If  the  crystal  is 
homogeneous  -  save  for  the  diffusing  species  -  this  assumption  amounts  to 
restricting  the  macroscopic  electric  field  to  values  too  small  to  effect  the 
jump  probability.  This  corresponds  to  the  assumption  that  ga  ^  «  W  where  a 
is  the  lattice  constant  and  W  is  the  barrier  height  if  the  jumps  are  thermally 
activated.  If  tunneling  predominates  the  situation  is  considerably  more  com¬ 
plicated  and  model  dependent.  Experimentally,  however,  this  appears  to  be  a 
good  approximation  at  least  at  high  temperatures.  With  this  simplification 
the  diffusion  equation  becomes: 

•  (7) 

C.  MACROSCOPIC  Field  Equations 

To  determine  4>(x,t)  we  impose  the  boundary  condition  that  <J>  »  0  at  the  surface 
and  make  use  of  Poisson's  equation  in  CCS  units. 

4irp 

£  ‘  (8) 

where  e  is  the  dielectric  constant  and  p  is  the  charge  density.  The 
interpretation  of  (8)  requires  some  care  in  that  the  potential  needed 
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in  (2)  and  (3)  is  the  potential  at  the  equilibrium  sites  of  the  diffusing  ion, 
and  hence  is  a  local  potential  which  may  be  subject  to  a  local  field  correction. 
If  so,  we  would  have: 

fLoc  =  %l  *  > 

where  is  the  field  at  the  ion  site,  ^  is  the  macroscopic  field,  £  is  the 
polarization,  and  a  is  the  appropriate  Lorentz  local  field  factor.  Now  in 
our  case  P  is  not  being  produced  by  the  diffusing  ions  themselves  and  hence  ? 
is  not  proportional  to  but  rather  to  E^.  The  proportionality  constant 

will  not  be  just  the  susceptibility  but  will  involve  Lorentz  local  field  factors 
appropriate  to  the  system.  The  net  result,  however,  is  that: 

?  • 

where  y  is  a  constant. 

Hence : 


(1  ♦  Y)E„ 


Thus : 


d<(> 


Loc 


dx 


“  *  «  TT 


and 

*Loc  =  (1  +  Y)  <J>M  ♦  6  . 

d4>. 

where  4.  is  macroscopic  potential  given  by  (8) ,  is  defined  by - =,  E 

LOC  '  qx 

and  6  is  a  constant.  The  requirements  of  charge  neutrality  for  the  crystal  as 
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a  whole  means  and  4>^oc  both  equal  0  at  the  surface.  Hence: 


Thus  <J>^  as  obtained  from  (8)  can  be  used  in  (7)  by  simply  replacing  q  in  (7) 
by  qeff.  However,  in  the  absence  of  experimental  evidence  to  the  contrary  we 


shall  suppose  that  y  ■  0  so  that  qe^  =  q,  and 


The  total  charge  density  is  made  up  of  tho  following  terms: 


PL  »  Charge  density  due  to  Lattice  Ions; 
p  =  -  n  e  »  Conduction  electrons; 

PT  =  -  n^e  =  "Trapped”  electrons; 

Pn  =  nq  =  Diffusing  ions; 

Pj  =  Non-diffusing  Ions; 

Py  a  Valence  electrons; 

where  -e  equals  the  charge  of  an  electron.  Noting  that  PL  +  Pv  =  0,  for  an 
insulator  we  have: 


p  *  pj  +  nq  -  e(nc  +  Rp) 


(9) 


In  the  absence  of  electric  fields, 


(10) 


where  N^  is  the  concentration  of  traps  with  trapping  level  E^,  8  =  1/kT,  and 
y  is  the  electron  Fermi  level.  Similarly: 
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I 


(11) 


n 

c 


E 


0(E) 


dE 


j^rr 7 


where  Ec  is  the  energy  of  the  bottom  of  the  conduction  band  and  D(E)dE  is  the 
number  of  states/cm^  between  E  and  E  ♦  dr. 

In  the  presence  of  a  macroscopic  electrostatic  potential,  equations  (10) 
and  (11)  take  the  form: 


(12) 


n‘(x)  ■  BtE)  TTgRbM  7  -  •  03, 

Ec  (x)  e  4  1 

In  practice,  nc(x)  is  always  small  enough  that  (13)  can  be  approximated  by: 
n  (x)  -  Nn  e-B[Ec (0)-u-e4>(x)]  0e*(x) 

U  S  '  (14) 

where  Ns  is  the  conduction  electron  density  at  the  surface  of  the  crystal  (where 
<Kx)  is  taken  to  be  zero). 

We  now  consider  a  special  case  of  particular  experimental  interest.  We 
suppose  that,  in  addition  to  the  conduction  band  (and  shallow  traps  which  can 
be  lumped  with  it),  we  have  deep  filled  traps  and  one  trapping  level  near  the 
Fermi  level.  This  case  can  be  made  to  fit  nearly  all  experimental  situations 
by  appropriate  choice  of  parameters.  The  concentration  of  electrons  on  the 
trap  near  the  Fermi  level  is: 
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"t  "  "'S(W  °-M 


e  ♦  1 


♦  1 


f  -Be*  } 

ie__  -  lj 

e 

**  J 

(15) 


where  is  the  density  of  the  trap,  Ej,  is  the  energy  of  the  trap  at  the 

surface  of  the  crystal,  and  N^.  is  the  concentration  of  electrons  on  the  trap. 
Hence : 


p(x)  ■  Px  ♦  n(x)q  -  cN$' e6e*(x)  -  e (NQ  ♦ 


(16) 


where  Np  is  the  concentration  of  filled  traps  and  Ns’  is  the  concentration  of 
electrons,  in  the  conduction  band  and  in  shallow  traps,  at  the  surface  of  the 
crystal.  Note  that  we  have  ignored  any  contribution  due  to  surface  states. 

These  can  be  included  by  an  additional  6  function  term  at  the  surface  if  needed. 

At  this  point  it  is  useful  to  note  some  of  the  special  cases  covered  by 
(16).  In  a  perfect  crystal  Pj  ■  NQ  -  N?  -  0  and  (16)  becomes: 

P(x)  ■  n(x)q  -  eNs'  e0®*^  .  (16a) 

In  the  case  of  several  deep  traps  well  below  the  Fermi  Level  (16)  becomes: 

P (x)  «  n(x)q  -  eNQ  -  eN.' e6®^30  .  (16b) 

Similarly,  most  experimental  situations  can  be  described  by  (16). 
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°-  T"°  Component  Diffusion  Equation  with  Macroscopic  Field 


He  are  now  In  a  position  to  write  down  a  system  of  equations  governing  the 
simultaneous  diffusion  of  H  and  D  in  the  presenee  of  a  macroscopic  electric 

field  they  help  to  produce.  We  assume  that  no  other  electrically  active  species 
is  diffusing.  From  (7)  and  (8)  wo  have: 


and 


3nH(x,t) 

92nH(x,t) 

— 5t — 

3nD(x,t) 

32nD(x,t) 

3t 

■g°-T7- 

j2<Kx,t) 

-4*P(x,t) 

dx2 

e 

given  by  equation  16. 

3 


kT  ST 


nH(x.t)  H £4 


nD(x.t)  ^ 


(17) 


(18) 


(19) 


- - -  AJIWCllUi 

of  the  sample  but  not  at  the  surface.  This  is  best  seen  by  considering  the 
li.if5y.no  in  (17)  and  (18)  above  (Pick’s  2nd  haw).  The  situation  of  interes 
experimentally  is  one  in  which  the  concentration  at  the  surface  is  constant. 
Hence  the  solution  of  Fick’s  2nd  Law  is  of  the  form 


n(x,t)  ■  Aq  ♦  I  An(t)  sin(n7rx/Jl)  , 
n*l 

where  i  is  the  thickness  of  the  sample.  Now,  strictly  speaking,  Fick's 
Second  Law  is  not  valid  at  the  surface  due  to  the  fact  that  the  crystal  ends 
and  hence  there  is  a  discontinuity  in  the  physical  parameters.  Mathematically 
we  add  an  "anti-crystal-  at  the  surface  and  extend  the  solution  to  negative  x 
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and  to  x  >  l.  Since  n(x,t)  -  A^  is  odd  at  x  *  0,  the  effect  is  to  make  both 
[n(x,t)  -  Aq]  and  d  n/dx  *  0  at  x  *  0,i.  Hence  the  discontinuity  can  either 
be  removed  by  adding  the  anti-crystal  or  can  simply  be  ignored  since  both  give 
precisely  the  same  answer  for  n  and  d2n/dx2  at  x  =  O.Jl.  In  other  words,  the 
symmetry  is  such  that  the  discontinuity  at  the  surfaces  can  be  neglected. 
However,  in  the  case  of  Eqs .  17-19,  this  is  no  longer  true,  (n  -  AQ)  and  <P  are 
odd  at  the  surface,  but  9/9x[n(d$/dx)]  is  not.  Thus  the  discontinuity  cannot 
be  ignored  and  we  must  simulate  the  surface  by  adding  the  term  required  to 
make  9n/9t  «  0  at  x  »  0,1, 

This  source  term  which  accounts  for  the  injection  (or  ejection)  of  ions 
at  the  surface  by  the  electric  field  is: 


as.  d 

‘  kT  37 


[<S(x)  4  6(x  -  Jt)] 


(20) 


where  6  is  the  Dirac  6  Function.  With  this  modification,  Eqs.  17-19  describe 
the  physical  system. 

We  now  need  appropriate  boundary  conditions.  We  shall  be  interested  in 
the  case  in  which  the  initial  concentrations  of  H  and  D  are  constant  throughout 
the  sample  except  at  the  surface,  and  the  surface  concentrations  are  constant  in 
time.  For  example,  we  will  have  conditions  of  the  type: 

njj  (x,0)  *  Nh  0  <  x  <  l 

nD  (x,0)  ■  Nd  0  <  x  <  l 

nH  (0't}  3  \  "  ”h  ^  (21) 

nD  (0»tJ  “  %  “  nD  (*’tJ 

<t>(0)  -  <|>(t)  -  0  . 
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\ 


Note  that  i-he  condition  on  <J>  is  simply  a  statement  of  overall  charge  neutrality 
for  the  crystal. 


E*  Small  Potential  Approximation 


The  system  Cl 7) -(19)  is  in  general  nonlinear  because  of  the  dependence 
of  p  on  <|>.  However,  if  0e<f>  «  1,  the  system  can  be  linearized.  The  resulting 
equations  admit  solutions  for  which  ge<f>  »  1.  Hence,  if  the  solution  does  indeed 
satisfy  the  condition  Be<p  «  1.  the  approximation  is  satisfactory.  In  this  spirit 
we  now  consider  p  in  this  limit.  Referring  to  Eq.  (16)  we  find: 


N„ 


*0 

e(ET-p) 


♦  i 


1  -  geifi 


“BCE--P) 
e  ♦! 


where 


%  0 
^  “B0L.-P) 
e  ♦  1 


1  ♦  ge<|> 


B(ET-p) 

6(et-p) 

e  ♦  1 


G  ♦  H<f>  , 


Hence 


where 


N_ 

_  *0 

i 

TO 

9 

o 

_ l 

B(et-p) 

© 

flter-V)  *  n  " 

e  ♦  1 

— 

e  ♦  1. 

p(x)  %  [px  -  eNs'  -  eND  -  ecj  ♦  nH(x)qH  ♦  nD(x)qD 

-  ♦[eH  ♦  Be2Ns'j  =  M  ♦  q^x)  ♦  q^fr)  -  N*(x)  , 

M  -|Pj  -  «  Ns»  ♦  Nd  ♦  gJ  ;  N  •  eH  ♦  0e2N5' 


(22) 


(23) 


(24) 

(24) 
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F. 


Solution  in  the  Small  Potential  Approximation 


To  solve  the  system  (17)- (19)  we  expand  the  concentrations  and  potential 
in  Fourier  Series  as  follows: 


00 


nHU,t)  *  Aq  ♦  Z  An(t)  sin(nwx/Jl) 
n»l 

1 

(25) 

00 

nI)(x,t)  =  Bq  ♦  Z  Bn(t)  sin(nirx/£) 
n=l 

1 

(26) 

00 

<j»(x)  ^  Z  Cn(t)  sin(nirx/£) 
n=l 

(27) 

Note  that  the  requirement  of  constant  n  and  <j>  at  the  surface  excludes  cosine 
terms  from  Eqs.  25-27.  Substituting  (25) -(27)  into  (19)  we  find: 


Hence : 


y  r  (W)2  ,  fmrx]  4tt 

'  c»  bn  slnm  ■  •  t 


00 

M  4  {,H*n  *  Vn  -  NC„>  *  Vo  *  V« 


Vo  *  Vo  ■  •  m  :  C„  (x)  ■  T  [Vn  *  Vn  -  NC„]  • 


Therefore 


,  WOOlA  ♦  qDB  ) 

n  (4n/e)N  ♦  (nir/4)2 


(28) 


(29) 


(30) 


Next,  substituting  (25)  and  (27)  into  (17)  we  find: 
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I  9A”  I’™']  .  „  r  i  MZ  j  fn**l  •  V|h  3 

TT  Slnm  -gH  An  W  S1"H  *  -W  57 


v 

nir 

.  fnirx) 

.  qHgH 

d 

w 

sinH 

kT 

3x 

oo 

z 


mir 


n,m*»l 


A *(*CJ  ‘-f  sin 


n '  m' 


nirx) 

3 

X 

00 

a  r  ^ 

mu 

fmirx) 

m cos 

m 

♦  Art  1  C 

0  ,  m 

m*l 

w 

cos(nrj 

r  A  fn7T]  cirifnlli) 

'  *H  MtJ  sln["T-J  TT~ 


r  7a  r  fm1r)  / fn1T)  fnirx)  fralTX]  fnn1)  .  fiMrx)  .  fi 

„ii  \a"c4tJ\W  cos  l-x-j  «sH  -  W  sinM  sinl' 


nnx 

nr 


fmir)2 


A°  Jl  ^ 


sin 


miTx 

nr 


(31) 


Note  that  (20)  insured  that  (dAQ/dt)  *  0.  Now: 


and 


cos(nr]  cos(nr]  "  j^os[(n  +  ®)  ♦  cos^(n  -  ®)  * 


sin 


fnirxl 

m 


sin 


mirx, 

nrl  B  2 


2  ^cos[(n  ■  *)  j]  '  cos  £(n  ♦  m)  . 


(32) 


(33) 


Substituting  (32)  and  (33)  into  (31),  multiplying  by  sin(mrx/A)  and  integrating 
over  x  from  0  to  A,  and  setting  qH  e  qD  «  e,  we  find: 


8An  .  fmt]2  epH  .  _  [mr] 

"St"  “  n  [t]  '  TT  VnlTJ 


fc]  -s-  z 

mA  C 

-n  p  *  *  r  P  -  n 

1 X J  im 

V  ;  P® 

p  m 

n  -  (p  ♦  a)  n  -  |p  -  m|2 

odd 

(34) 
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Similarly, 


SB 

T5T  -  -  *db„ 


2  e«L 

TJ  ‘  TT  B0Cn 


nir 

T 


28r 


mri  e  _ 

■T*  r—  £  rnA  C 

■J  ™  pm  P  m 
odd 


1  n 


5-EJ-!L-7  * 

-  (P  ♦  m)  n  -  |p  -  m | 


(35) 


The  system  (30),  and  (34) -(35)  cannot  be  solved  analytically  but  it  is  very 
easily  done  numerically.  Thus  as  soon  as  experimental  numbers  (work  of  Sec.  Ill) 
are  available  for  the  parameters  the  problem  will  be  completed. 


V.  Determination  of  Conduction  Band  Parameters 

Considerable  experimental  work  has  been  done  on  applying  the  theory  of 
hydrogen  solubility  of  the  last  progress  report  to  the  determination  of  con¬ 
duction  band  parameters,  trapping  levels,  etc.  The  difficulties  encountered 
in  this  work  demonstrate  quite  compellingly  precisely  why  these  measurements 
are  so  essential.  The  first  parameters  which  must  be  determined  are  the  proton 
binding  energy  and  the  conduction  band  density  of  states.  After  these  have  been 
determined  to  sufficient  accuracy,  it  will  then  be  possible  to  measure  the  Ti 
interstitial  binding  energy  and  to  start  measuring  thermal  trapping  energies 
for  electrons  on  various  defects  and  impurities,  and  to  calculate  the 
equilibrium  concentration  of  both  <H>  and  <Ti>  (H  and  Ti  interstitials). 
Determination  of  <H>  binding  energy  and  conduction  band  density  of  states 
requires  measurement  of  equilibrium  concentration  of  <H>  and  conduction 

electrons,  over  a  fairly  wide  temperature  range,  with  P„  and  P  held 

H,  H-0 
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constant,  in  a  concentration  ange  where  the  electron  Fermi  level  is  high  enough 
that  we  can  be  sure  that  all  electron  traps  arc  filled.  In  order  for  our  re¬ 
sults  to  be  meaningful,  however,  we  must  simultaneously  hold  the  concentration 
of  <Ti>  at  a  very  low  level,  since  we  have  no  quantitative  means  of  measuring 
their  concentration.  Fortunately,  the  theory  is  sufficiently  detailed  to 
provide  a  variety  of  checks  to  determine  whether  or  not  this  requirement  ' 
being  met.  So  far,  we  have  encountered  considerable  difficulty  in  selecting 
operating  conditions  which  yield  adequate  <H>  concentrations  without  significant 
<Ti>  concentrations.  This  same  difficulty  has  been  the  source  of  many  incorrect 
interpretations  in  published  work.  The  most  recent  measurements  indicate  that 
we  now  have  the  problem  under  control,  however,  and  quantitative  results  should 
be  obtained  in  the  next  few  weeks. 

Because  of  the  extreme  importance  of  these  measurements  (severa’  other 
projects  can  be  done  properly  only  after  these  measurements  are  completed), 
most  cf  our  efforts  recently  have  been  devoted  to  this  probiw .  If  we  are 
successful  in  conq>leting  these  determinations,  it  will  be  the  first  time, 
to  our  knowledge,  the  defect  and  impurity  thermodynamics  of  a  system  such  as 
rutile  have  been  worked  out  theoretically  and  the  necessary  experimental 
parameters  obtained.  This  will  make  possible  a  wide  variety  of  experiments 
that  could  not  previously  have  been  done,  and  will  be  a  major  step  forward 

in  obtaining  a  useful  characterization  of  this  system. 
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VI.  Far  IR  Spectrometer 


The  spectrometer  is  essentially  complete  and  the  first  experiments  are 
now  being  planned.  We  vrill  begin  with  a  general  survey  of  transmission  and 
reflection  properties  with  special  attention  being  devoted  to  looking  for  the 
pla* ma  frequency. 
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